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Starting from a covariant formalism of the Sunyaev-Zeldovich effect for the thermal and non- 
- - - thermal distributions, we derive the frequency redistribution function identical to Wright's method 

assuming the smallness of the photon energy (in the Thomson limit). We also derive the redistri- 
bution function in the covariant formalism in the Thomson limit. We show that two redistribution 
functions are mathematically equivalent in the Thomson limit which is fully valid for the cosmic 
microwave background photon energies. We will also extend the formalism to the kinematical 
Sunyaev-Zeldovich effect. With the present formalism we will clarify the situation for the discrep- 
' ancy existed in the higher order terms of the kinematical Sunyaev-Zeldovich effect. 
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I. INTRODUCTION 



o 
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The Sunyaev-Zeldovich (SZ) effect [U, i, i, Q, which arises from the Compton scattering of the cosmic microwave 
Q ■ background (CMB) photons by hot electrons in clusters of galaxies (CG), provides a useful method for studies of 
cosmology. For the reviews, for example, see BirkinshawH and Carlstrom, Holder and ReeseQ. The original SZ 
^2 , formula has been derived from the Kompaneets equation(7[ in the non-relativistic approximation. However, recent 
^ X-ray observations (for example, Schmidt et al.[8] and Allen et al.Q) have revealed the existence of high-temperature 
CG such as fc^Te ~20keV. Wright 10] and Rephaeli and his collaborator pH [T^ have done pioneering work including 
' the relativistic corrections to the SZ effect for the CG. 

In the last ten years remarkable progress has been made in theoretical studies of the relativistic corrections to 
the SZ effects for the CG. Stebbins[3] generalized the Kompaneets equation. Challinor and Lasenby 14] and Itoh, 
Kohyama and Nozawajl^] have adopted a relativistically covariant formalism to describe the Compton scattering 
^vq ■ process and have obtained higher-order relativistic corrections to the thermal SZ effect in the form of the Fokker- 
Planck approximation. Nozawa, Itoh and Kohyama; 16] have extended their method to the case where the CG is moving 
with a peculiar velocity with respect to the CMB and have obtained the relativistic corrections to the kinematical 
SZ effect. Their results were confirmed by Challinor and Lasenby [13] and also by Sazonov and Sunvaevfisl [l9|. Itoh, 
] Nozawa and Kohyamafg^ have also applied the covariant formalism to the polarization SZ effect [1,[3]- Itoh and his 
collaborators (including the present authors) have done extensive studies on the SZ effects, which include the double 
scattering effect (2l|, the effect of the motion of the observer 1221 ■ hig h precision analytic fitting formulae to the direct 
numerical integrations [23l[23] and high precision calculations [2 5l [26| . The importance of the relativistic corrections is 
^ ' also exemplified through the possibility of directly measuring the cluster temperature using purely the SZ effect [27]. 

On the other hand, the SZ effect in the CG has been studied also for the non-thermal distributions by several 
groups [2^ [29l. [30| . The non-thermal distribution functions, for example, the power-law distributions, have a long tail 
in high electron energy regions. Therefore the relativistic corrections for the SZ effect could be more important than 
the thermal distribution. 

Shimon and Rephaeli 31] have discussed on the equivalence of different formalisms to the SZ effect. The relativistic 
SZ effect has been studied analytically so far in three different approaches. The first method is the calculation of the fre- 
quency redistribution function in the electron rest frame using the scattering probability derived by Chandrasekhar [s^ . 
This method was used by Wright [l3] and extended by Rephaeli [ll[. We call it as Wright's method in the present 
paper. The second approach solves the photon transfer equation in the electron rest frame. This approach was 
used by Sazonov and Sunvaevp^. We call it the radiative transfer method. The third approach is the relativistic 
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generalization of the Kompaneets equationp?], where the relativistically covariant Boltzmann coUisional equation is 
solved for the photon distribution function. This approach was used by Challinor and Lasenbv|il4j and Itoh, Kohyama 
and Nozawaflq. We call it as the covariant formalism in the present paper. In Shimon and Rephaelifsij they have 
shown the equivalence between Wright's method and the radiative transfer method. They also have claimed the 
equivalence between Wright's method and the covariant formalism. However, no mathematical relations are shown 
between the redistribution function in Wright's method and the expression of the scattering probability in the covari- 
ant formalism. Therefore their claim is incomplete. In the present paper we will show explicitly that two approaches 
are mathematically equivalent. 

On the other hand, recently Boehm and J. Lavalle[30| also have discussed the equivalence of the different approaches 
for the SZ effect in the non-thermal distribution. They have shown that the radiative transfer method is equivalent 
to the covariant formalism. However, they have concluded that Wright's method is incorrect. In the present paper 
we will show that their conclusion is incorrect. We will show that Wright's method, which has been widely used in 
the literature, is still fully valid. 

The fourth method for the study of the SZ effect is the direct numerical integration of the rate equation of the 
photon spectral distortion function. The first-order calculation in terms of the optical depth r was done by Itoh, 
Kohyama and Nozawap^ for r <IC 1. The full-order calculation was done by Dolgov et al.^] for t 1. The rate 
equation in the present formalism has a simple form. Therefore it is more suitable for the direct numerical application. 
We will present the numerical calculation elsewhere [3^. 

The present paper is organized as follows. In § II, we show the equivalence between Wright's method and the 
covariant formalism of the SZ effect for both thermal and non-thermal distributions. We also derive the rate equations 
and their formal solutions for the photon distribution function and for the spectral intensity function. In § HI we extend 
the formalism to the kinematical SZ effect, and derive the rate equations in Wright's method. Finally, concluding 
remarks are given in § IV. 



II. SUNYAEV-ZELDOVICH EFFECT 



A. Equivalence between Covariant Formalism and Wright's Method 

Let us consider that both of the CG and the observer are fixed to the CMB frame. As a reference system, we choose 
the system which is fixed to the CMB. (Three frames are identical in the present case.) In the CMB frame, the time 
evolution of the photon distribution function n{Lu) is written as follows |15l|: 

Mi^) '^^P d^p'd'^k'W {n{ij)[l + n{io')]fiE)~n{ij')[l + 72{Lu)]f{E')} , (1) 



dt J (27r 

2ujLo'EE' ' ^ ' 

^ = -(4 + -)+4^^(- + ^)'-4-^(- + A), (3) 

K^-2{p-k) ^ -2ljE{1- ^^i) , (4) 
k' ^2{p-k') ^ 2uj'E{l-Pfi') , (5) 

where e is the electric charge, m is the electron rest mass, W is the transition probability of the Compton scattering, 
and f{E) is the electron distribution function. The four-momenta of the initial electron and photon are p = (E,p) 
and k = (oj, k), respectively. The four-momenta of the final electron and photon are p' = {E' ,p') and k' = (oj', k'), 
respectively. In Eqs. ^ and (O, P = IpI/E, h = cos9 is the cosine between p and fc, and fi' — cos9' is the cosine 
between p and k' . Throughout this paper, we use the natural unit h = c = 1, unless otherwise stated explicitly. For 
later convenience we rewrite Eq. ([3]) as follows: 

X^Xa + Xb, (6) 

(7) 

(k-k'f 

Xb = -4^ . 8 
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By eliminating the (5-function, Eq. ([T]) is rewritten as follows: 



X + n(c^')K(i^) - n{u')[l + n(c^)K(i^')} , (9) 

dr = TieOTdt , (10) 

7 = , (11 

f{E)^n,7r^p,{E), (12) 

where rig is the electron number density, ctt is the Thomson scattering cross section, and Pe{E) is normalized by 
dpp^pe{E) — 1. By choosing the direction of the initial electron momentum (p) along z-axis, the photon momenta 

k and k' are expressed by 

fc = (^Vl - /"^cos(/)fc, v^l - M^sin^fe, /x^ , (13) 

k' = u' (^y/l- /i'2cos0fc' , Vl - Ai'^sin^fc' , ^i') , (14) 

where (pk and t/ifc' are the azimuthal angles of k and k', respectively. Inserting Eqs. (I13p and (fT4|l into Eqs. ^ and 
dl]), one obtains 

^ _ (1 - cose)^ 1 - cose 

^74(l-/?;,)2(i-/3^')' 72(1-/3m)(1-/3;x')' ^ ^ 
V ^/'^V^' (1-cose)^ 

^ ^ (1-/3m)(1-/3/^')' ^ ' 



1 — + (ti>/7m)(l — COS0) ' 

C0S8 = + Vl - ^1 _ cos{(f)k - 4>k' ), (18) 

where cos0 is the cosine between k and k' . It should be noted that Xa and will not be mixed each other under 
an arbitrary Lorentz transformation, because Xa depends only on /i, /3 and 7, whereas Xb depends also on oj and 
uo'. 

Now let us introduce the transformations for /i and /i' which will play a key role in the present paper. 

-Mo + /3 i.^. 

1 - |3^l'Q 

where /iq = cos6'o and /ig = cos^q are cosines in the electron rest frame. The suffix denotes the electron rest 
frame throughout this paper, unless otherwise stated explicitly. Equations ((T9)) and (pO)) are the composition of the 
Lorentz transformation for the photon angles from the CMB frame to the electron rest frame and the transformation 
9q vr — 6*0, 6'q TT — Oq. Note that the latter transformation is not essential. Applying Eqs. (fT9| and (|20p to 
Eqs. - (HH]), one obtains as follows: 

Xa = 1 + cos'^Qq, (21) 



X ^f^y^' (l-cosGo)^ 
^ {^mj UJ {1 - (3fio){l ^ Pfi'o) ' ^ ' 

^ ^ 1 -/^Mo (-23) 

Lu 1 — (3^0 + {uj/^m){\ — cos8o) 

cosGo = hoh'q + sJ I - iil \J I - fi'^ cos{(j)k -ipk'), (24) 
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where cosOq is the cosine between k and k' in the electron rest frame. It can be seen that Eq. ((2T|) was surprisingly 
simplified compared with Eq. p3|) . On the other hand Eq. ((^ did not change its form compared with Eq. (fTB]) . As 
will see later in this section, Eqs. (|2ip and ([22|) are the key points for connecting the covariant formalism with Wright's 
method. The terms Xa and Xb did not mix each other by the above reason. Furthermore Xa is the expression in 
the electron rest frame, whereas Xb is not, because it contains /? and jm. 
The phase space volumes are transformed as follows: 

^ 72(i_'^^^)2 ^"^'o > (26) 

where (Ppo = p^dpd^od(j)p, dflk'o = d^Qd(l)k'- Note that the z-axis was chosen along k direction for the d^p integration. 
With these variables Eq. ^ is re-expressed by 



dn{uj) 



d'po / dQk'o^- — [-] X 



2 



dT 647r2 7 " 74l-/3/io(l-/9/io)' 

X {n{Lo)[l + n{u;')]p,{E) - n{Lu')[l + n{u;)]p,iE')} . (27) 

In deriving Eq. ((27)) we used the relation 7^(1 — Pfi) = (1 — Pfio)^^- 

Before proceed to the next step, some explanations might be necessary for Eq. ^7\ . In Eq. ([?7)) photon zenith 
angles (/io and /Iq) are described in the electron rest frame with the transformations of Eq. (fT9|) and ([20]) . On the 
other hand, energies (w, to' and p) and azimuthal angles {(j^k, (j^k' and (/)p) are left in the CMB frame. As seen later in 
this section, this peculiar hybrid coordinate system makes the connection from the covariant formalism to Wright's 
method in a straightforward manner. It is needless to say that the familiar Klein-Nishina formula in the electron rest 
frame will be obtained by the Lorentz transformations to = luqj (1 — /3/io) and lu' = ujq"/ (1 — /S/ig) and inserting into 
Eqs. (EI]) and (122D. 

Now let us introduce an assumption which was also used in Boehm and Lavalle[30|. 

7— < 1 . (28) 
m 

For the CMB {kBTcMB = 2.348 x IQ-^'eV) photons w < 5 x IQ-^eV is well satisfied. Then Lo/m < 1 x 10"^ 
which implies 7 ^ 10^. Therefore as far as the CMB photon energies are concerned, Eq. (pS)) is fully valid from 
the non-relativistic region to the extreme-relativistic region for the electron energies. With Eq. (|28p the following 
approximations are valid. 



UJ 1 - f3fio ' 
Xb = 



UJ 

7— 



(29) 
(30) 

O (/37-)l , (31) 

(E') = (E) [ 0[TcMB/Te)] for thermal distribution , . 

Pe\ ) Pe\ [1 + 0(70;/™)] for power law distribution ' ^ ' 

As seen from Eqs. (pn)) ~ ([5^ . the Thomson limit is realized in the scattering kinematics by the assumption of Eq. 
With these approximations Eq. (|27|) is reduced to 



E' = E 



^^6^1 / ^"^'V (l-L)^ (^ + - ■ (^^^ 

Furthermore the cjjk' -integral can be performed and one obtains 

d^k'{l + cos^eo) = l + M^ + l{l-^^lKl-t^^)■ (34) 

Inserting Eq. ([M]) into Eq. ([55]) and assuming the spherical symmetry for pe{E), one obtains as follows: 



dn{uj) 
~d7~ 



dpp Pe{E) 



1 rl 1 



^ d^o y ^ ^^0 2^4 (1 _^^q)3 -^(^"' ^0) ["(^^0 " "(^)] > (35) 
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/(Aio,Mo) = 

According to Wright [To*] we introduce a new variable s by 



^ u/ ^ 1-M ^ 
u 1 - /3/io ' 

which imphes dfiQ = — (1//3)(1 — f3fio)e^ds. Then Eq. ([35|) is finaUy rewritten by 



9t 



P(s,/3) 



rtJ.2{s) 

4 / c?Mo(l - /3mo) 



1 



(1 - Pf^oY 



f/(MO,Mo) 



where 



W -ln[(l+/3)/(l-/3)] 

Mf,-[l-e^(l-/3Mo)]//3, 
-1 



[l-e-^(l+/3)]//3 



for s < 
for s > 

for s < 
for s > 



(36) 

(37) 

(38) 
(39) 

(40) 
(41) 

(42) 
(43) 



Equation (j39p is the probabihty for a single scattering of a photon of a frequency shift s by an electron with a velocity 
/3, which is described in the electron rest frame. By using the identity relation 1 — /3^o = g'*!! ~ Pf^o), Eq- (EH) is 
identical to P{s; (3) (Eq. (7)) in Wright [lOj. Thus Wright's redistribution function has been derived from the covariant 
formalism. 

Now we will derive the redistribution function in the covariant formalism under the assumption of Eq. (j28p (the 
Thomson limit). The derivation is straightforward but lengthy. We will give the derivation in Appendix A and will 
quote the result here. The expressions which correspond to Eqs. psp and ([39]) in the covariant formalism (in the 
CMB frame) are 



dn{Ld) 
dr 



2/372 



dsP{s,(3) [n(e"w) - n(w)] , 



dfi'f{^,fi') , 



/(m,a^') 



'Ml(s) 

[1-6^(1-/3//)]//?, 



- 2 



1 - nn' 



72(1-/3m)(1-/3m') 



(44) 

(45) 
(46) 

(47) 



where Smax, /^i(s) and /i2(s) are defined in Eqs. (|40|) . (|42|) and (|43|) . respectively. In the present paragraph we show 
that P{s,(3) is identical to P(s,/3). In order to show the equivalence, we apply the transformations of Eqs. (fTOll and 
(HOI) to Eq. (gni). First, inserting Eqs. ^ and ^ into Eq. (gT]), one obtains 



/(m,m') = .f(A*o,/^o) • 



The variables /i' and fiQ have the relation 



which implies 



dn' 



7^(1-/3^0)^ 



(48) 



(49) 



(50) 
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and boudary values 



_ I fi2is) at fi' = ^i(s) , . 

^ /.i(s) at m' = M2(s) ■ ^^^^ 



Inserting Eqs. (|i5|) - ([5T|) into Eq. (^5]) . one finally obtains 



gS /•M2(s) 1 



which is identical to Eq. ([55)1 . Therefore one obtains 

P(s,/3)=P(s,/?). (53) 

Thus the equivalence between the covariant formalism of the Boltzmann collisional equation fl5| and Wright's 
method p^. [ll[ has been shown mathematically under an assumption ^oj/m ^ 1, where the assumption is fully 
valid for the CMB photon energies. It should be emphasized that no non-relativistic approximations are made for the 
electron energies in deriving Eqs. ([55)1 and ([IS]) . This is the reason why the calculations by two different formalisms 
produced same results for the SZ effect even in the relativistic electron energies. In Appendix B, we have also shown 
the derivation of Eq. ((27l) in terms of the Klein-Nishina cross section formula. 

Boehm and Lavalle.30j also discussed the equivalence between the radiative transfer approach and the covariant 
formalism. However, they concluded that Wright's method was incorrect. We conclude that their conclusion is 
incorrect. The reason why they lead the erroneous conclusion is as follows. They start with the covariant form for 
the squared Compton amplitude (their Eq. (43)). They derived the familiar Chandrasekhar's form (their Eq. (50)) by 
taking the non-relativistic limit (/3 0) in their Eq. (49). Because of the non-relativistic approximation they used, 
they concluded that Wright's method (Eq. (50)) should not be used for the relativistic calculation. On the other 
hand, we have also started with the same covariant form for the squared Compton amplitude. We have derived the 
same expression (Eq. (34)) without taking the non-relativistic limit. We have shown that Eq. (34) is connected to 
its covariant form by the Lorentz transformations of Eqs. (19) and (20). Therefore Wright's method is equivalent 
to the covariant formalism. We conclude that their criticism is incorrect. Shimon and Rephaehim also claimed the 
equivalence between the covariant formalism and Wright's method. Their Eq. (19) looks similar to Eq. ([55]) . however, 
no mathematical relations are shown explicitly in their paper between W in their Eq. (19) and P{s; (3) of Wright . 

B. Rate Equations and Formal Solutions 

We now proceed to derive the rate equations and their formal solutions. Since two formalisms are equivalent, one 
can use either P(s,/3) or P(s,/3). We start with Eq. ([55)1 and rewrite as follows: 



dn{ijj) 



dsPi{s)[n{e''uj) ^n{Lo)] , (54) 



1 

2„,5 ~ 



dr 

P,{s)^ I dp0'j'p,{f3)P{s,p), (55) 
/3™„ = (l-e-l^l)/(l + e-l^l), (56) 



where Pe{f3) = m^pe{E). As seen from Eq. ([55|l . Pi{s) is the probability for a single scattering of a photon of a 
frequency shift s averaged over the electron distribution function, which is so called the redistribution function of a 
shift s. The total probability is dsPi{s) = 1. Multiplying oj^ to Eq. ([54)) . one obtains the rate equation for the 
spectral intensity function. 



dljiu) 
dr 



dsPiis) [e"3"/(e"w) - I{uj)] , (57) 



where I(lu) — uj^n{uj)/2Tr'^ is the spectral intensity function for uj. Now let us introduce the following key identity 
relations: 

P(s,/3)e-3^ = P(-s,/3), Pi{s)e-^' = Pi{-s) . (58) 
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The derivation is straightforward. Inserting Eq. (|58p in Eq. (|57|) and replacing s by —s, one obtains the rate equation 
for the spectral intensity function. 



dr 



dsPi{s) [l{e-'u;)-I{Lj)] . (59) 



It should be remarked that 71(6" uj) appears in RHS of Eq. whereas /(e ^w) appears in RHS of Eq. ([55)1 . It is 

also straightforward to show that Eq. (I54p satisfies the photon number conservation. 



d 

d7 



dujuj^niuj) = . (60) 



Let us now derive formal solutions for the rate equations Eq. ((54|) and Eq. ([59|) . We consider an ideal condition 
that the CG is infinitely large. We introduce a new function n{LU,T) by 

?i(a;) = e^^ri(cj, r) . (61) 

By inserting Eq. (|6ip into Eq. ((54| . one obtains the equation for n(uj, t). 

dfi{ijj, t) 



- , dsPiis)n{e'Lu,T), (62) 
where dsPi{s) = 1 was used. Equation (|62p can be integrated and one has 

fi{uj,T) = no{Lj) + dX dsPi{s)il{e''uj,X) . (63) 



In deriving Eq. (j63p an initial condition 71(0;, r = 0) = no(w) was used, where no(w) is the initial photon distribution 
function. We solve Eq. (|63p with a successive approximation method. The first-order term is obtained by inserting 
no(w) into RHS of Eq. ((63)) . 



hi{u} , t) = no{Lu) + T / c?sPi(s)no(e''w) . (64) 
The second-order term is also obtained by inserting hi{uj,T) into RHS of Eq. ((63l) . 

/oc 
(isFi(s)no(e'*cj) 
-00 

+ ^ / rf^^-P2(s)no(e^o;), (65) 

/OO 
dsiFi(si)Pi(s-si), (66) 
-00 

where P2{s) is the probability (redistribution function) of a shift s for the double scattering. By repeating the above 
procedure iV -|- 1 times, one obtains the (TV -t- l)-th order term. 



nN+i{^,T) = no{uj) + 2_^— / dsP,-(s)no(e''[j) , (67) 
j=\ ' •'-00 

/OO /"OO 
dsiPiisi)--- ds,_iPi(s,_i)Pi(s-^s,), (68) 
-00 J —OO 

where Pj{s) is the probability (redistribution function) of a shift s for the multiple scattering of the j-th order. By 
taking the limit A*" — > cxd in Eq. ((67)) and replacing limjv^oo '^Ar(^, t) = n{uj, r), one finally obtains the formal solution 
for n{uj). 

/OO 
dsP{s,T)no{e''uj) , (69) 
-OG 
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Multiplying uj^ to Eq. (|69|) and using P{s,T)e = P{—s,t), and also replacing s by — s, one obtains the formal 
solution for I{uj). 

/oo 
dsPis,T)Ioie-'Lu), (71) 
-oo 

where Iq{uj) = [x)'^/(27r^)no(w). Note that this solution can be also derived directly from Eq. ([55)1 . Note also 
that Eq. (j70p is the Poisson distribution function. The distribution function is commonly used, for example, in 
Birkinshaw[5]. In the present paper, however, Eq. ()70p is derived as a natural consequence of the present formalism. 

In practical cases, the CG has a finite size and the optical depth is small (r <^ 1), therefore the first order 
approximation is sufficiently accurate for the study of the SZ effect. From Eqs. (I5n | - (|7T|) one obtains the following 
familiar forms for the distortion functions. 

An(uj) = n(uj) — no(w) 

dsPi{s)[noie'Lo)-no{co)] , (72) 

A/(t^) = /(w)-/o(^) 
« r 



dsPiis) [loie-'u) - loiuj)] , (73) 
r = <TT J dirie . (74) 
The integral in Eq. (f74|) is done over the photon path length in the CG. 



III. KINEMATICAL SUNYAEV-ZELDOVICH EFFECT 

Let us now consider the case that the CG is moving with a peculiar velocity (3c {—Vc/c) with respect to the CMB. 
As a reference system, we choose the system which is fixed to the CMB. The z-axis is fixed to a line connecting the 
observer and the center of mass of the CG. (We assume that the observer is fixed to the CMB frame.) In the present 
paper we choose the positive direction of the z-axis as the conventional one, i.e. the direction of the propagation of 
a photon from the observer to the cluster, which is opposite to that of Nozawa, Itoh and KohyamaJBj. In the CMB 
frame, the time evolution of the photon distribution function n(uj) is same as for the thermal SZ effect as shown in 
Nozawa, Itoh and Kohyamaflij. They are given by Eqs. HI)-®. The electron distribution functions are Lorentz 
invariant and are related as follows: 

f{E) = UEc) , (75) 
J{E')^UE',), {76) 



Ec - S7c [l + Pc-f3) , (77) 
E'c = S'tc (l + Pc- /3 ') , (78) 



7c = , , (79) 

where the suffix c denotes the CG frame. Therefore the formalism of § II will be directly applicable to the present 
case. A modification should be made to the electron distribution fmiction Pe{E) by 



Pe [E) = pe,c [E-fc \l+(3c-(3\) , (80) 

where Pe,c{Ec) is normalized by dpcPcPe,c{Ec) = 1- To proceed the calculation, one expresses the product (3c ■ (3 
in the coordinate system where k is parallel to the z-axis. Then one obtains 

(3c- 0= (3c(3 {A*c/i + ^/l - - fi^ cos{(()c - 0p)} , (81) 
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where /ic and 0c are the cosine of the zenith angle and the azimuthal angle of /3c, respectively. By applying the 
transformation of Eq. to Eq. (|8ip . one obtains 



Pep 



1 - f3fj.o 



; cos 



(82) 



Inserting Eqs. (|80p and 
Wright's method. 



into Eq. (|33p . one obtains the expression for the CG with non-zero peculiar velocity in 

dn{uj) 



dr 



dpp / d/io 



1 



1 



/(Aio,/^o) 



■27r 



d(t)pPe.c (^E^c 



(83) 



Shimon and Rephaeliim also obtained the expression for the kinematical SZ effect based upon Wright's method, 
which is similar to Eq. For the expression of /3c ■ P, Eq. (|82p agrees with their Eq. (39). As discussed in their 



paper, however, they have an extra factor 7c ^1 + /3c • f3j in Eq. ([55]) which comes from E^/E in their phase space 

factor, see their Eq. (37). As discussed also in Nozawa, Itoh, Suda and Ohhata[26j. the reason of the discrepancy 
is because they used the phase space in the CG frame instead of the CMB frame. As far as the present formalism 
is concerned, we have used the CMB frame as a reference system. Therefore there are no extra factors needed in 
Eq. 183]) . We conclude that the result of Shimon and Rephaeli is in error by the extra factor. 

Let us now proceed with Eq. For most of the CG, /3c -C 1 is realized. For example, /3c ~ 1/300 for a typical 

value of the peculiar velocity Wc^lOOO km/s. In Nozawa, Itoh and Kohyama[lB| they made an expansion in terms of 
/3c in the Fokker-Planck approximation. They found that O{0^) terms are negligible for most of the CG. Therefore 
we will keep 0{(3c) terms and neglect higher-order terms in the present paper. In this approximation the electron 
distribution function is approximated as follows: 



PeAEc) ~ PeiE) { 



1 - 



a/3c • /3 
E 



knT, 



for Pf, [E) cx p " 

for pe{E) oc E-" 

for Pe{E) (X exp^—E/ksTe 



(84) 



For simplicity, we consider the thermal distribution function. (Only a minor modification will be needed for the 
power-law distributions.) Inserting Eq. (|82p into Eq. (j84p the integral for the azimuthal angle is performed. 



1 

2^ 



2tt 



E-fc 



1 + /3c • /3 



^^PeiE) 



(85) 



where 9^ = ksTe/m. Repeating the same procedure done in § II, one obtains the rate equations for the case of the 
CG with nonzero peculiar velocity. 



dn{uj) 

dr 
dljuj) 

dr 



dsPi (s, l3c,z) [nie^oj) - n{oj)] , 

DO 
C 

Pi{s, I3c,z) ^ Pi{s) + Pc,zPi,k{s) , 



(86) 
(87) 



where Pi{s) is Eq. ((55)) and /3c. z = PclJ'c is the peculiar velocity parallel to the observer, because the photon direction 
is along z-axis. In Eq. ([SS]) . Pi^k{s) is the redistribution function due to the peculiar velocity of the CG. It is given 

as 



Pi,k{s)^ [ dP(3^j^p4P)PK{s,P)., 



1 



M2(s) 

(1-/3mo)^ 



/(a*o,Mo) . 



(89) 
(90) 
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where /Xq, /^i(s), fJ.2{s) and Pmin are defined in Eqs. (|4T|) . (|42|l. (|43)) and (|56|) . respectively. It should be remarked 
that Eq. ([57]) is expressed by e~^''/(e''w) instead of /(e^'^tj) in Eq. ([51]). This is because P(s,/3)e~^'* = P{-s,/3) as 
shown in Eq. ([58|l . however, Pk{s, (3)e~^'' ^ Pk{—s,(3). For the power-law distributions, {j/de) should be replaced 
by a//3^ and a in Eq. (|90p for the p-powcr distribution and the i?-power distribution, respectively. 

Finally, one obtains the distortions of the photon spectrum and the spectral intensity in the first order approxima- 
tion. 

/oo 
dsPi {s, f3c,z) [noie'uj) - 7io(c^)] , (91) 
-OO 

/OO 
dsPi{s,(3,.,) [e-'^'Ioie'uj) - /qH] . (92) 
-OO 



IV. CONCLUDING REMARKS 

We started with a covariant Boltzmann coUisional equation of the SZ effect shown in Itoh, Kohyama and Nozawa[l^ 
for thermal and non-thermal distributions. First we have applied a rational transformation (Eqs. (fT9|) and (pO| ) to the 
photon angles, which is essentially a Lorentz transformation for photon angles from the CMB frame to the electron 
rest frame. The transformation has made the expression for the transition probability surprisingly concise form. Then 
we have introduced an assumption used by Boehm and Lavalle[30|. namely ■yiv/ni <C 1 (the Thomson limit). The 
assumption is fully valid for the CMB photon energies. Under the assumption, we have derived the redistribution 
function P(s,/3), which is the probability for a single scattering of a photon of a frequency shift s by a electron with 
a velocity /3. The obtained redistribution function is identical to that of derived with Wright's methodfiol. [Tl|. 

Similarly, starting from the covariant Boltzmann collisional equation of the SZ effect for thermal and non-thermal 
distributions, we have derived the redistribution function P{s, j3) in the covariant formalism under the assumption 
^Lo/m <^ 1. We have shown that P{s,P) is identical to P(s,/3). They are connected by the Lorentz transformation 
of Eqs. (fTO|) and (|^D|) . Thus we have shown mathematically that Wright's method is equivalent to the covariant 
formalism under the assumption 'yto/m <^ 1. This result guarantees that existing works which used Wright's method, 
for example, BirkinshawiSj, Enfilin and Kaiser [2^ and Colafrancesco et al.[2^, are still fully valid. This result also 
explains the reason why two different calculations for the thermal SZ effect agree extremely well even for the relativistic 
electron energies. 

We have also extended the present formalism to the kinematical SZ effect. Starting from the covariant Boltzmann 
collisional equation for the kinematical SZ effect, we have repeated the same procedure. We have derived the redistri- 
bution function for the CG with nonzero peculiar velocity in Wright's method. We have compared the present result 

with that of Shimon and Rephaeli[3l|. The obtained redistribution function is differ by a factor 7c (^1 + Pc ■ fi^ ■ We 

have clarified the discrepancy between their result and others[lB, [13, [3- Their result is in error by the factor. 
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APPENDIX A: REDISTRIBUTION FUNCTION IN COVARIANT FORMALISM 

In this appendix we will derive the redistribution function in the covariant formalism. The starting equation is 
Eq. dS]). 

dn{uj) 3 /■ ,3„ /■ ,o 1 1 /^^"^ 



X {n{u;)[l + n{uj')]p,iE) - n{uj')[l + n{uj)]peiE')} . (Al) 
Then we assume the Thomson limit ^Lu/m <^ 1, which implies the approximations 

^ « 1^ (A2) 
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and Xb ^ 1, £" « i? and Pe{E') « Pe{E). Under the assumption, Eq. (|A1|) is approximated as 



dr 647r2 7 ^^<=v ^ ^2 (1 - /3^')2 

where is given by Eq. (jlSp . In Eq. (jA3p the 0^/ -integration can be done as 



1 



27r 



XAdcpk' = 2 



(i-/Vr + -(i-M')(i-M'') 

1 - M/^' 



72(1 -/3m) (1-/3^0 

Assuming the spherical symmetry for pe{E), Eq. (jASP is further simphfied. 



dn{uj) 
dr 

/(m,m') 



dpp Pe(£') / / 



-1 J-i 



\ , 1 I -[3(1 



272(1-/3^0 



7Ti/(m,m')["(^') -"(t^)], 



(1-;.m')'+2(1-m')(1-m") 
74(l-/3/i)'(l-/3M')2 



1 - M/^' 



72(1-/3m)(1-/3m') 



Now let us introduce a new variable s by 



uj' 1 - /3a< 



1 - /3At' ' 

which implies dfi = — (1//3)(1 — j3ii')e^ds. Then Eq. (jA5[) is finally rewritten by 

dn{uj) 



dr 



dpp^Pe{E) / dsP{s,l3)[n{e''uj)-n{uj)] 



2s rP-2{s) 



W J,, is) 

[l-e^(l-/3M')]//3 



dn'f{n,fi') , 



(A3) 



(A4) 

(A5) 
(A6) 



(A7) 

(A8) 

(A9) 
(AlO) 



where Smax, fJ-iis) and ^2(5) are given by Eqs. (|40| . (|42| and ((43|) . respectively. Equation (|A9|1 is the redistribution 
function in the covariant formalism, which is described in the CMB frame. 



APPENDIX B: KLEIN-NISHINA CROSS SECTION 

In this appendix we will derive Eq. (|27p in terms of familiar Klein-Nishina cross section formula. Notations are 
same as those in the main text, unless otherwise stated explicitly. As a reference frame we choose the electron rest 
frame. The energy-momentum conservation gives the relation for the photon energies as follows: 



1 



u!q 1 -I- (wo/m)(l - cos6o) ' 



(Bl) 
(B2) 



cosGo = MoMo + \/l - Mo \/l - M? cos(0fco - (t)k'„) , 
where Qq is the scattering angle. The Klein-Nishina cross section formula in the electron rest frame is expressed by 



da 



1 



sin" 80 



where is the classical electron radius. With Eq. (|Bip one obtains the following usefuU relation. 



(B3) 



(B4) 
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Inserting Eq. (jB4p into Eq. (jB3p one can rewrite the Klein- Nishina formula as follows: 

2 



da 



1 



l + cos^Go+f^^) i^:^(l-coseo) 



(B5) 



It is needless to say that one obtains the Thomson cross section by taking the limit too/m <^ 1 and uj'q/luq —^ 1 in 
Eq. (lB5l) . 

Now let us introduce the transformation from the electron rest frame to the CMB frame, where the electron is 
moving with a velocity /3. The photon energies cu and uj' in the CMB frame are related to and uj'q by the Lorentz 
transformation 



uj = ujol (1 - ^3^io) , 
c.'=47(l-M) , 

where fiQ = cosSq and /Iq = cos^q. With the variables oj and oj' one obtains 

2 / /\ 2 " 



da 
dU,u 



1 



1 - /?Mo 



1 + cos^ Oo 



\ ^ oj' (1 — cos 0o)^ 
7m y w (1 - /3^o)(l - /3Mo) 



(B6) 
(B7) 



(B8) 



As seen from Eq. (jBSj) the square bracket in the RHS is identical to Xa + ^s, where they are defined by Eqs. (PT|) 
and (211). Note that Eq. (jB8|) is the expression in the hybrid coordinate system, where the energies are described in 
the CMB system, whereas the zenith angles are described in the electron rest frame. 

The cross section is defined by the transition rate divided by the flux of the incident particles. The flux in the CMB 
frame is 



p ■ k 



(B9) 



Therefore, one can write Eq. Ilj in terms of the cross section in the CMB frame as follows: 

dn{uj) 



dt 



= -2 



-^(1-/3.)^ — 



{n{u)[l + n{u;')]}{E) - n{uj')[l + n{u)]f{E')} . 



(BIO) 



Since the cross section is Lorentz invariant, one can rewrite Eq. (|B10|) with the Klein-Nishina cross section in the 
hybrid system of Eq. (|B8p as follows: 



dn{uj) 
dt 



d^k 



(27r)3' '"^'ydnk'^ 

{n(c^)[l + n{^')]f{E) - + n{uo)]f{E')} 

Rewriting the phase space volume d^p by 

1 



d-'p 



7'(1-/3mo)2 

and inserting Eqs. (jBSp and (jB12p . one finally obtains 



(Bll) 



(B12) 



dn{uj) 
dr 



647r2 



74l-/?/io (1-M)2 



X {n{u)[\ + n{uj')]pe{E) - n{uj')[l + n{uj)]p,{E')} 



(B13) 



In deriving Eq. (jB13p we used the relations 7^(1 — = (1 — /9/io) f{E) ~ ir'^nePeiE), dr = Ueardt and 
aT = S-K/irl- One finds that Eq. ((BT3)) is identical to Eq. (P7)) . 
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